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The Figure of Merit Science Working Group (FoMSWG) recently forecast the constraints on dark 
energy that wiU be achieved prior to the Joint Dark Energy Mission (JDEM) by ground-based 
experiments that exploit baryon acoustic oscillations, type la supernovae, and weak gravitational 
lensing. We show that cluster counts from on-going and near-future surveys should provide robust, 
complementary dark energy constraints. In particular, we find that optimally combined optical and 
Sunyaev-Zel'dovich effect cluster surveys should improve the Dark Energy Task Force (DETF) figure 
of merit for pre- JDEM projects by a factor of two even without prior knowledge of the nuisance 
parameters in the cluster mass-observable relation. Comparable improvements are achieved in the 
forecast precision of parameters specifying the principal component description of the dark energy 
equation of state parameter as well as in the growth index 7. These results indicate that cluster 
counts can play an important complementary role in constraining dark energy and modified gravity 
even if the associated systematic errors are not strongly controlled. 



I. INTRODUCTION 

Counts of galaxy clusters are a potentially very power- 
ful technique to probe dark energy and the accelerating 
universe (e.g. [2 H d H [51 [51 H [i i]). The idea is an 
old one: count clusters as a function of redshift (and, 
potentially, mass) , and compare to theoretical prediction 
which can be obtained either analytically or numerically. 
Recently, Rozo et al. [in] have obtained very interesting 
constraints on erg from Sloan Digital Sky Survey (SDSS) 
cluster samples using the relation between mass and opti- 
cal richness (the number of red-sequence galaxies in the 
cluster above a luminosity threshold). This follows re- 
cent dark energy constraints using optical [TT] and X-ray 
observations of clusters [T2j [131 HI] • 

In this paper we calculate the potential of cluster 
counts to improve combined constraints from the other 
three major probes of dark energy: baryon acoustic oscil- 
lations (BAO) , type la supernovae (SNIa) and weak grav- 
itational lensing (WL) . We are motivated by the recently 
released report of the Figure of Merit Science Working 
Group (FoMSWG; (15j) that studied and recommended 
parametrizations and statistics best suited to addressing 
the power of cosmological probes to measure properties 
of dark energy. While the FoMSWG report was mainly 
aimed at figures of merit to be used in the upcoming 
competition for the Joint Dark Energy Mission (JDEM) 
space telescope, the applicability of its results and rec- 
ommendations is general. 

We address quantitatively how ongoing and upcoming 
cluster surveys, in particular the South Pole Telescope 
(SPT, [IS]) and the Dark Energy Survey (DES, [TT]), 
can strengthen the combined "pre- JDEM" constraints on 
dark energy considered in the FoMSWG report — that 
is, combined constraints expected around the year 2016. 
To model cluster counts, we utilize recent results from 
Cunha |18j which optimally combine future optical and 
Sunyaev-Zel'dovich (SZ) observations of clusters to esti- 



mate the constraints on dark energy. 



II. INFORMATION FROM CLUSTER COUNTS 
AND CLUSTERING 

The subject of deriving cosmological constraints from 
cluster number counts and clustering of clusters has been 
treated extensively in the literature (see e.g. [T51 [151 
[m mi). In this work we use cross-calibration for 
two observable proxies for mass: Sunyaev-Zel'dovich flux 
(henceforth SZ), and optical observations — which iden- 
tify clusters via their galaxy members — (henceforth 
OPT). While we focus on optical and SZ surveys, our 
results are applicable to combinations of any cluster de- 
tection techniques. In particular, planned X-ray surveys 
such as eRosita |13], WFXT gl], and IXO gS] will have 
mass sensitivity competitive with, and complementary 
to, the SZ and optical surveys. 

Our approach closely follows that in [TH] and we refer 
the reader to that publication for basic details. In brief, 
cluster counts in a bin of the observables are calculated 
by integrating the mass function dn/dM over mass, vol- 
ume, and the observable proxy in the appropriate range. 
We adopt the Jenkins mass function in this work, though 
results are weakly dependent on this choice. Clustering 
is given by the sample covariance of the mean counts 
in different redshift bins. The contribution of clustering 
to the constraints is very small when cross-calibration 
is used [15]. We allow for scatter in both the relation 
between mass and the observable proxy, and the rela- 
tion between true and estimated photometric redshifts. 
Results from both simulations (e.g. [261 [22]) and observa- 
tions (e.g. [lHlllSl[3n]) suggest that the mass-observable 
relations can be parametrized in simple forms with log- 
normal scatter of the mass-observable about the mean 
relation. Other works (see e.g. [31]) suggest that the dis- 
tribution of galaxies in halos may be more complicated. 
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We assume lognormal scatter for the mass-observable re- 
lation as well as for the photometric redshift errors. We 
have neglected any theoretical uncertainties in the mass 
function, galaxy bias, and photometric redshifts, all of 
which must be independently known to a few percent so 
as not to affect cosmological constraints [21} [32] . 

We fix the photo-z scatter to ~ 0.02, the expected 
overall scatter of cluster photo-z's in the Dark Energy 
Survey |17j . Our "theorist's observable" quantity, which 
we feed into the Fisher matrix formalism to obtain con- 
straints on cosmological parameters, is the covariance of 
the counts — defined as the sample covariance plus the 
shot noise variance — in different redshift bins. 

We adopt the same surveys and parametrizations de- 
scribed in [18], namely, an SZ and an OPT survey on the 
same 4000 sq. deg. patch of sky. 

Let Afobs be the observable proxy for mass (from ei- 
ther SZ or OPT survey) . For the SZ survey, we define the 
Afobs threshold for detection to be Af*^ = lO^^' ^/i^i ATq, 
complete up to z = 2, based on the projected sensitivity 
of the South Pole Telescope. We parametrize the mass 
bias (the difference between the true mass and SZ Afobs), 
and the variance in the mass-observable relation, respec- 
tively as 



X^M^'^'^z) = lnAfo'"''"-|-ailn(l-|-z) (1) 
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Fiducial mass values of all nuisance parameters are zero, 
except for the scatter which is set to ctq — 0.25 in the 
fiducial model. Our choice of scatter is somewhat con- 
servative given recent studies which suggest that (Tq < 0-2 
(see e.g. |33jV However, hydrodynamics simulations by 
[M] find a scatter of (-1-32%, —16%) about the median for 
clusters of M ^ 3.0 x 10^^ Af©, which matches our choice. 
In total, there are six nuisance parameters for the mass 
bias and scatter (InAfp"'''; ^i; '^'h^ ^i)- 

For the optical survey the mass threshold of the ob- 
servable is set to Af'h = IO^^-^H^^Mq and the redshift 
limit is z = 1, corresponding to the projected sensitiv- 
ity of the Dark Energy Survey. Different studies suggest 
a wide range of scatter for optical observables, ranging 
from a constant a\nM = 0.5 [22] to a mass-dependent 
scatter in the range 0.75 < crinM < 1-2 [3S]. Using weak 
lensing and X-ray analysis of MaxBCG selected optical 
clusters, Ref. ^ estimated a lognormal scatter of ~ 0.45 
for P(Af|A^obs), where M was determined using weak 
lensing and Afobs was an optical richness estimate. We 
choose a fiducial mass scatter of ai^M = 0.5 and allow 
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FIG. 1: First three principal components for the OPT-I-SZ 
cluster survey with Planck priors. Solid lines refer to the case 
with no prior on the 16 nuisance parameters, while the dashed 
lines correspond to the case of perfectly known nuisance pa- 
rameters. 

for a cubic evolution in redshift and mass: 

InA/i^'"'^ (Afobs, ^) = InA^'"^ + ai ln(l + z) 

-I- a2(lnA'fobs - InAfpivot) (3) 

3 

t^ln M (^^obs ,Z) = O-g -t- ^ 6, Z' 

1=1 

3 

+ Ec'(lnA//obs-lnAfpivot)' (4) 

i=l 

We set Afpivot = lO^^/i'^Af©. In all, we have 10 nuisance 
parameters for the optical mass errors (InAfQ'^*^, ai, 02, 
Cq, 6i, Ci). 

There are few, if any, constraints on the number of pa- 
rameters necessary to realistically describe the evolution 
of the variance and bias with mass. Ref. [20j shows that 
a cubic evolution of the mass-scatter with redshift cap- 
tures most of the residual uncertainty when the redshift 
evolution is completely free (as assumed in the Dark En- 
ergy Task Force (DETF) report [37]). Note too that we 
employ more nuisance parameters to describe the optical 
survey than the SZ survey because the former is expected 
to have a more complicated selection function. For the 
cross-calibration analysis, we assume the correlation co- 
efficient between optical and SZ scatter p, defined in [18] , 
to be fixed to zero; the same paper shows that the cross- 
calibration results are insensitive to the value of p for 
PG hl,0.6]. 

In total, we use 6-1-10=16 nuisance parameters to de- 
scribe the systematics of the combined OPT+SZ cluster 
survey. While generous, this parametrization assumes 
a lognormal distribution of the mass-observable relation 
that may fail for low-masses. We have also implicitly 
assumed that selection effects can be described by the 
bias and scatter of the mass-observable relation. By the 
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year 2016, we expect significant progress in simulations 
of cluster surveys that will allow us to better parametrize 
the cluster selection errors. 



III. COMPLEMENTARY PROBES AND 
FIGURES OF MERIT 

To model the power of complementary probes of dark 
energy, we adopt the pre-JDEM information (that is, 
combined information projected around year 2016) based 
on estimates of the Figure of Merit Science Working 
Group [131. These estimates include information from 
BAG, SNIa, WL and the Planck GMB sateUite. We use 
these probes in combination, without or with clusters. 
Note that systematic errors have been included in all of 
these methods (see Ref. [15]). 

The FoMSWG figures of merit are described in the 
FoMSWG paper [TS] and we review them here very 
briefly. There are a total of 45 cosmological parame- 
ters, 36 of which describe the equation of state w{z) 
while the others are mostly standard cosmological pa- 
rameters (plus a couple of nuisance ones that have not 
been explicitly marginalized over). Gne flgure-of- merit 
is the area in the WQ-Wa plane [371 131], where w{z) — 
wq + Wa{^ — a) — Wp + Wa{ap — a) and where Wp and 
flp are the "pivot" parameter and the scale factor; we 
adopt FoM = l/{(T{wp) x <7[wa))- The growth of den- 
sity perturbations is described by a single parameter, the 
growth index 7, which is a free parameter in the fltting 
function for the linear growth of perturbations [39]. The 
flgure-of-merit in the growth index is simply its inverse 
marginalized error, 7F0M = 1/(7(7). 

A much richer (and less prone to biases) description of 
the equation of state is achieved through computing the 
principal components (PCs) of dark energy [ID], ei{z) 

35 

1 w(a) = ^ aiej(a), (5) 

i=0 

where Ui are coefficients, and ei(a) are the eigenvectors 
(see [TS] for details). The associated flgure of merit con- 
sists of presenting the shapes ei{z) in redshift and com- 
puting the associated accuracies a{ai) with which the 
coefficients can be measured [TS] . 

Combining the different cosmological probes is 
achieved by adding their associated Fisher matrices. We 
add the 45 x 45 Fisher matrix for clusters (marginal- 
ized over the mass nuisance parameters) to the combined 
BAG+SNIa-l-WL+Planck Fisher matrix and report the 
improvement in the flgures of merit and accuracies in the 
PCs as well as shapes of the new PCs. 



IV. RESULTS 

Gur baseline is the combined pre-JDEM 
BAG-l-SNIa-l-WL+Planck case from the FoMSWG 



report. The baseline uncertainties in various dark 
energy parameters, after marginalization over all nui- 
sance and cosmological parameters, are cr(wo) = 0.10, 
a{wa) = 0.31, a{wp) = 0.028 {zp = 0.40), FoM = 116, 
and cr(7) = 0.21 (7F0M = 4.8). These constraints 
are dominated by BAG+Planck which alone yield 
cr(wo) = 0.15, a{wa) = 0.44, (j{wp) = 0.037, FoM = 61. 
In comparison, WL-|-Planck and SNIa-|-Planck yield 
FoM = 9.8 and 0.42, respectively. 

We flrst consider the cluster information alone, with 
only a Planck prior adopted from [15 . The constraints 
in this case, assuming flat external priors on cluster nui- 
sance parameters, are cr(wo) = 0.10, a{wa) = 0.41, 
a{wp) = 0.036 {zp = 0.28), FoM = 66, and (7(7) = 0.17 
(7F0M — 6.0). These constraints are comparable (and 
complementary) to the BAG-|-Planck constraints. Fig. [l] 
shows the first three principal components for the com- 
bined (GPT-I-SZ) cluster survey combined with Planck. 
Two cases are shown: completely unknown and perfectly 
known nuisance parameters. In the first case, the first 
principal component peaks at z ~ 0.3, which is not sur- 
prising given that most clusters are at z < 1: while this 
peak sensitivity is at lower redshifts than that for BAG 
surveys, it is at slightly larger z than the peak for SNIa. 

As seen in Fig. [T] adding priors to nuisance parame- 
ters moves the cluster PC weights to higher z. This is 
easy to understand: freedom in the nuisance parameters 
has progressively more deleterious effects as redshift in- 
creases, as can be deduced from Eqs. Q-Q. Thus, pri- 
ors on the nuisance parameters restore the ability of the 
survey to probe higher redshifts, and push the principal 
components to higher z. 

Next, we combine the cluster information with 
BAG-l-SNIa-l-WL-l-Planck. The amount of informa- 
tion that clusters contribute is a strong function of 
the systematics assumed, in our case, a function of 
the external priors on the nuisance parameters. We 
find that clusters provide very significant improve- 
ment in the figures of merit even with uninformative 
(flat) priors on the cluster nuisance parameters. The 
new clusters-l-SNIa-l-WL-l-BAG-l-Planck flgure of merit 
is 206, which is nearly a factor of two better than 
BAG+SNIa+WL+Planck alone. The pivot error is 
a{wp) = 0.022 {zp 0.41) which is ~ 25% better. Con- 
straints on Wq and Wa improve by about 50% each, since 
with clusters cr(wo) = 0.065 and <j{wa) = 0.214. Con- 
straints on growth improve by more than a factor of 2, 
with (7(7) = 0.099 (7F0M = 10). The main effect of 
including clusters is to provide signiflcant additional in- 
formation on the growth of density perturbations. Gf 
the complementary techniques we consider, only weak 
Icnsing probes the growth, but our flducial cluster model 
gives (slightly) stronger constraints on growth than the 
pre-JDEM combination of BAG-t-SNIa-|-WL-|-Planck. 

The left panel of Fig. [2] shows the best-determined 
three principal components for the flducial pre-JDEM 
survey adopted from |15| when clusters are added, i.e. 
clusters-l-SNIa-l-WL-l-BAG-l-Planck. For the cluster sur- 
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FIG. 2: Left panel: First three principal components for the SNIa+BAO+WL+Planck pre-JDEM combination alone (solid 
curves) and for the same combination with the addition of clusters (dashed curves). For the latter we assume OPT+SZ cluster 
survey with flat (i.e. uninformative) external priors on nuisance parameters. Right panel: Uncertainty in eigencoefflcients of 
w(a), (T(ai), for the SNIa+BAO+WL+Planck pre-JDEM combination alone (circular points), and for the same combination 
with the addition of clusters (square points). We also show the ratio of the improvement in each eigencoefBcient when clusters 
are added (dashed line - scaled down by a factor of 10 for clarity). 




0.001 0.01 



0.1 1 10 100 

on mass scatter 



o 



o 









OPT 




SZ \ 




OPT+SZ 


1 





0.001 0.01 



0.1 1 10 

on mass scatter 



100 



FIG. 3: Improvement in the figures of merit relative to the pre-JDEM combination (BAO+SNIa-|-WL-|-Planck) when cluster 
information is added. In both panels we add a uniformly increasing prior on the nuisance parameters as explained in the 
text, and show (on the x-axis) the effective resulting prior on the scatter in mass CTinM(Af, z), that is, the (square root of the) 
left-hand sides of Eqs. ([2]| or Q that correspond to the prior values of the nuisance parameters on the right-hand sides at 
z — 1.0 and M — \{i^^h~^ Mgj. For the OPT-I-SZ case, the uncertainty in the optical scatter is shown on the x-axis. The 
left panel shows the improvements in the (DETF) FoM, while the right panel shows improvements in the inverse error in the 
growth index, 7F0M. 



vey we assume OPT+SZ with flat external priors on nui- 
sance parameters. We see that clusters make the total 
principal components look more "cluster-like" (compare 
to Fig. [1]) since they add a lot of information to the to- 
tal. For the same reason, clusters move the weight of the 
best-determined PC toward lower redshifts. 

The right panel of Fig. [2] shows the accuracies a(a.i) 
with which the coefficients on of the principal components 
can be measured. The contribution of clusters becomes 
more pronounced for higher PC number, leading to a 
nearly constant fractional improvement in error. With 
clusters, the fourth eigencoefhcient is about as well con- 



strained as the second eigencoefficient in the baseline case 
without clusters. 

If external priors on the nuisance parameters are avail- 
able, the full power of cluster constraints is even more 
evident. To model such priors, we scale all errors by the 
same fractional value — for each nuisance parameter pi, 
we let Fa — > Fii( \ -\- a) where a varies from zero (flat 
prior) to inflnity (sharp prior). Consequently, the addi- 
tional information in each nuisance parameter is a fixed 
fraction of the original (unmarginalized) error in the pa- 
rameter. While other choices are possible for adding pri- 
ors, we settle on this simple prescription to illustrate the 
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effects of external information on nuisance parameters. 

The left panel of Fig. |3]shows the ratio of the figure of 
merit that includes all probes, FoM, and the FoM with 
clusters left out (that is, FoM/FoMq where FoMq = 116). 
We consider three cluster survey scenarios: an optical 
survey, an SZ survey with optical follow-up for photomet- 
ric redshift measurements only ^ , and the cross-calibrated 
OPT-I-SZ survey. On the x-axis we show the effective 
prior on the scatter in mass (the quantity {crf^j^f{z))^^'^ 
from Eqs. ([2| and ^) — priors are added to other nui- 
sance parameters as well, we simply do not show them. 
The plot shows that the OPT-I-SZ combination improves 
the total FoM by more than a factor of four if the scat- 
ter is known to high precision. In the more realistic cases 
where mass scatter (and other corresponding parameters) 
are known to finite accuracy from independent measure- 
ments, we still see improvement by factors of ^ 3. 

Priors on nuisance parameters contribute to the infor- 
mation content only if they are substantially stronger 
than the intrinsic ( "self-calibrated" ) uncertainties in 
these parameters; for the scatter in mass, for example, 
this impHes the knowledge of (o'f^Mi^))^^^ better than 
0(1) as Fig. [3] shows. Comparing the OPT, SZ, and 
OPT-I-SZ cases, we see that as prior information ap- 
proaches zero (high values of CTprior), the cross-calibration 
provides a lot of extra information relative to OPT or SZ 
alone. 

The right panel of Fig. [3] shows the corresponding 
figure of merit for the growth index 7. Even stronger 
improvements are now seen, with the 7 figure of merit 
increasing between a factor of two (flat priors) and ten 
(infinitely sharp priors). However, we caution that simu- 
lations of modified gravity models need to be done to de- 
termine whether the impact of modified structure growth 
on the cluster abundance is adequately captured by the 
7 parameter. Nevertheless, the right panel of Fig. [3] 
indicates that clusters appear to have at least as much 
potential to improve the pre-JDEM constraints on the 
growth history of the universe as they do for the expan- 
sion history (a similar conclusion has been reached in 
Ref. [41] for a specific modified gravity model). We also 
see that the SZ survey is more useful for improving 7F0M 
than the DETF FoM; this is because SZ probes higher 
redshifts, which allows for improved constraints on the 
redshift evolution of the growth of structure and hence 
7- 



V. DISCUSSION: IMPLICATIONS OF THE 
ASSUMPTIONS 

In this section we discuss the validity of the assump- 
tions we made and the consequences of varying those 



^ For SPT, optical follow-up is expected from the DES, the Blanco 
Cosmology Survey (BCS) and the Magellan Telescope. 



assumptions. We divide our assumptions into optimistic 
and pessimistic. 

The assumptions we consider optimistic are: 

• The optical mass threshold (M*'' = IO^^-^/i^iMq); 

• The SZ mass threshold (M*!^ = 10^^-^ h'^MQ); 

• Perfect selection for both SZ and optical cluster 
finding; 

• SPT area (4000 sq. deg.; could be less); 

• Known functional form of the scatter in the 
mass-observable relation (lognormal); 

• No mass dependence in the SZ mass-observable 
scatter (see Eq. (|2|). 

• Perfect knowledge of photometric redshift errors. 

The assumptions that are arguably pessimistic are: 

• No other cluster techniques (e.g. X-ray or weak 
lensing) are available to further cross-calibrate 
cluster counts; 

• Large fiducial value of scatter for both optical 
(cto = 0.5) and SZ (ctq = 0.25); 

• Area of DES (4,000 sq. deg.; could be as large as 
10,000 sq. degrees); 

• Low redshift range of optical cluster-finding 
(z<l); 

• Cubic polynomial evolution of redshift scatter 
for optical and SZ and mass evolution of optical 
scatter (Lima & Hu, 2005 [20] show that cubic red- 
shift evolution of the scatter yields near-maximal 
degradation of cosmological parameters) ; 

• Constraints are based on our current knowledge of 
cluster physics, while the field is developing rapidly. 

The first three optimistic assumptions are the most 
important. Since the mass- function falls rapidly with in- 
creasing mass, the lower mass bins contain most of the 
clusters, and are therefore most relevant. Cunha (2009) 
[351 shows that cross-calibration decreases the sensitivity 
of the constraints to the mass threshold somewhat. Here 
we have checked that increasing the optical limit from 
logM"^ = 13.5 to 13.7, or the SZ hmit from 14.2 to 14.5 
degrades the figures of merit by 10-20%. Increasing both 
leads to 30% degradation in the FoMs. The importance 
of uncertainty in photometric redshift errors has been 
studied extensively by [21]. For the surveys we consider 
here, it is not unreasonable to assume that large enough 
training sets will be available to sufficiently constrain the 
evolution of the redshift errors and characterize the sur- 
vey selection. 

Of the pessimistic assumptions, the first one is espe- 
cially significant: for example, if X-ray information is 
available (as expected from surveys such as eRosita), 
then X-ray plus optical cross-calibration alone can lead 
to excellent dark energy constraints even with the unex- 
pected failure of one or more of our SZ assumptions. 

To test assumptions about the functional form of the 
scatter, we added another Gaussian to both optical and 
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SZ mass-observable relations^, for a total of 27 new pa- 
rameters (43 total) ; the new parameters describe the evo- 
lution with redshift and mass of the mean and variance 
of the new Gaussian (cf. Eqs.[3]and|4|, the ratio between 
the two Gaussians describing each mass-observable rela- 
tion, and the correlation coefficient between optical and 
SZ (see [55^). The figures of merit degrade by merely 15- 
20%. The small additional degradation is a consequence 
of the fact that the new nuisance parameters do not in- 
troduce significant new degeneracies with cosmological 
parameters. If we instead add 4 parameters to charac- 
terize the mass dependence of the SZ scatter and bias, 
the degradations are even weaker, being < 5%. Intu- 
itively, adding mass-dependent evolution of the SZ bias 
and scatter is not as important as the functional form of 
the scatter because the SZ probes too narrow a range of 
masses for the evolution to be significant. 

The arguments and tests outlined in this section show 
that the assumptions made in this paper are not overly 
optimistic, and that the unforeseen systematic effects 
would have to be rather capricious in order to lead to 
significant further degradations in the cosmological con- 
straints. 



VI. CONCLUSIONS 

We have shown that galaxy clusters are a potentially 
powerful complement to other probes of dark energy. As- 
suming optimally combined optical and SZ cluster sur- 
veys based on fiducial DES and SPT expectations and 
allowing for a generous set of systematic errors (a to- 
tal of 16 nuisance parameters), we have shown that the 
constraints on the figure of merit expected in 2016 from 
baryon acoustic oscillations, type la supernovae, weak 
lensing, and Planck improves by nearly a factor of two 
when clusters are added. This improvement is achieved 



without any external prior knowledge on the cluster 
mass-observable nuisance parameters (but also without 
explicitly allowing for errors in the theoretically predicted 
mass function or cluster selection). 

We have further illustrated the cluster contribution to 
constraints by computing, for the first time, the principal 
components of the equation of state of dark energy for 
clusters alone and clusters combined with other probes. 
We found that the first cluster principal component peaks 
at z ~ 0.3, indicating the "sweet spot" of cluster sensi- 
tivity to dark energy. This redshift increases slightly if 
external information on the cluster nuisance parameters 
is available. Each eigencoefficient of the principal compo- 
nent expansion is improved by about a factor of two when 
clusters are added, indicating that the improvements ex- 
tend to well beyond one or two parameters. 

Finally, we have shown that measurements of the 
growth index of linear perturbations 7 (which is a proxy 
for testing modified gravity) improve by a factor of sev- 
eral with cluster information. While this particular calcu- 
lation depends on assumptions about the modified grav- 
ity model, it broadly illustrates the intrinsic power of 
clusters to measure growth and distance separately and 
to obtain useful constraints on modified gravity explana- 
tions for the accelerating universe. 

We conclude that cross-calibrated cluster counts have 
enough intrinsic information to significantly improve con- 
straints on dark energy even if the associated systematics 
are not precisely known. 
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